National Aeronautics and 
Space Administration 

Dryden Flight Research Center 

P.O. Box 273 

Edwards, California 93523-0273 



Reply to 

Attn of: Dryden Technical Publications April 1, 2004 


TO: All Holders of NASA/TP-2003-212034, Dated December 2003 

FROM: T/Technical Publications Office 

SUBJECT: Errata Sheet for NASA/TP-2003-212034, Dated December 2003 


NASA/TP-2003-212034, Aging Theories for Establishing Safe Life Spans of Airborne 
Critical Structural Components, by William L. Ko has an error in table 2 on page 14. Please 
make the following changes to this document. 

1. Delete the current text on page 14. 

2. Insert the new attached page 14 with the corrected column C in table 2. 

3. Staple the errata sheet to the inside front cover of the document. 

Thank you for your cooperation. 


Michael H. Gorn 
Acting Chief of Code T 


Attachment 
as stated 


(45) 


«W)‘ 


1 — — 


(7 

V max/ 


( 2) 1 ” 
C(m - 2) 


[ AM k 




In equation (45), the right-hand side is known, and the left-hand side has two unknowns: <7 max and o s . Thus, 
there are multiple combinations of c max and o s for finding the equivalent-constant-amplitude load spectrum. The 
most practical way is to choose the value of the mean stress a s of the random loading spectrum by inspection, and 
then calculate the equivalent maximum stress G from equation (45). 


NUMERICAL EXAMPLE 


The numerical example is to demonstrate how to calculate the number of safe flights for a flight test program. 
The example chosen is the NASA B-52B carrier aircraft pylon hooks (one front hook, two rear hooks) (figs.l, 2) 
carrying the Pegasus winged three-stage solid rocket (44,629 lb) up to high altitude (approximately 40,000 ft) for 
air-launching, firing and sending the payload into orbit. 

Material Properties 

The material properties of the B-52B pylon front and real - hooks are listed in table 2. 


Table 2. Material properties of B-52B pylon hooks. 


Part name 

Material 

°U 

ksi 

<*Y 

ksi 

T U 

ksi 

Kr 

L c 

ksi /im 

C 

in. ,, ■ r . — 
(ksiVm.) 

cycle 

m 

n 

Front hook 

Inconel 718 

175 

145 

135 

125 

9.220 X 10“ 12 

3.60 

2.16 

Rear hooks 

AMAX MP35N 1 " 

250 

235 

141 

124 

2.944 X 10 -11 

3.24 

1.69 


* Inco Alloys International, Inc., Huntington, West Virginia, 
t H. C. Starck, Inc., Cleveland, Ohio. 


Input Crack Data 

The fictitious surface crack at the critical stress point of each hook is assumed to be semicircular in shape 
( c = a, fig. 14). This is based on the observations of the past failed B-52B pylon two rear hooks. The crack 
propagation initiation sites of the failed hooks were almost semi-circular surface cracks of depths a = 0.03 1 in. and 
a = 0.038 in. respectively for the inboard and outboard hooks (ref. 2). Note from figure 14 that the semicircular 
surface crack has the relatively high value of Q. Other crack parameters used were: A = 1.12 for the surface crack, 
and M K = 1.0 for the high (hook depth) to (crack depth) ratios. 
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ABSTRACT 


New aging theories have been developed to establish the safe life span of airborne critical structural 
components such as B-52B aircraft pylon hooks for carrying air-launch drop-test vehicles. The new aging theories 
use the “equivalent-constant-amplitude loading spectrum” to represent the actual random loading spectrum with 
the same damaging effect. 

The crack growth due to random loading cycling of the first flight is calculated using the half-cycle theory, and 
then extrapolated to all the crack growths of the subsequent flights. 

The predictions of the new aging theories (finite difference aging theory and closed-form aging theory) are 
compared with the classical flight-test life theory and the previously developed Ko first- and Ko second-order 
aging theories. The new aging theories predict the number of safe flights as considerably lower than that predicted 
by the classical aging theory, and slightly lower than those predicted by the Ko first- and Ko second-order aging 
theories due to the inclusion of all the higher order terms. 

NOMENCLATURE 


A crack location parameter 

a depth of semi-elliptic surface crack, in. 

a° c operational limit crack size, in. 

a P initial fictitious crack size established by proof load test, in. 

a l crack size at the end of the /-th flight, in. 

1 n / — — m 

C material constant in Walker crack growth rate equation, — (ksi Vim ) 

cycle 

c half length of surface crack, in. 

E complete elliptic function of the second kind 

Fj number of safe flights at the first flight 

F i number of safe flights based on Ko first-order aging theory 

F i number of safe flights based on Ko second-order aging theory 

* 

F number of safe flights based on Ko discrete aging theory 

* 

F| number of safe flights based on Ko closed-form aging theory 

/ fraction of proof load (f < 1) 

ksi kilopounds 

K i mode I critical stress intensity factor, ksi Vim 

K m mode I stress intensity factor associated with a , ksi Vin. 

AK mode I stress intensity amplitude, ksi Vim 

i, j = 1, 2, 3, ... , integers associated with the i- th, or j- th half cycle 

k modulus of elliptic function 



/ 


= 1, 2, 3, , integer associated with the /-th flight 

M, flaw magnification factor 

m Walker exponent associated with maximum stress o , or stress amplitude ( cr max - t7 j ) 

N, number of stress cycles consumed during the /-th flight 

^ total total number of stress cycles allowed for safe flights 

n Walker exponent associated with stress ratio R 

Q surface flaw and plasticity factor 

R stress ratio, /? = C7 min /C7 max 

SRB/DTV solid rocket booster/drop test vehicle 

V A front hook load, lb 

V bl left rear hook load, lb 

V B r right rear hook load, lb 

V P A front hook proof load, lb 

V P BL left rear - hook proof load, lb 

V P BR right rear - hook proof load, lb 

Act i amount of crack growth induced by the /-th flight, in. 

8a- crack growth induced by the /-th half cycle, in. 

6 c angular location of critical stress point measured from horizontal axis, degree 

(Tj tangential stress at critical stress point of front hook, ksi 

G-, tangential stress at critical stress point of left rear hook, ksi 

o\ tangential stress at critical stress point of right rear hook, ksi 

* * * * 

a = (T| or o 2 or a 3 , ksi 

❖ p 

(i| value of (T| associated with V A , ksi 

❖ P 

g 2 value of g 2 associated with V BL , ksi 

❖ p 

g 2 value of c 3 associated with V BR , ksi 

Gjj ultimate tensile stress, ksi 

Gy yield stress, ksi 

°max maximum stress, ksi 

G m | n minimum stress, ksi 

G s mean stress, ksi 

(T ; tangential stress at hook inner boundary, ksi 

(ad maximum value of G . , ksi 

t max 1 
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ultimate shear stress, ksi 


% 

( ) ■ quantity associated with the i-th stress half cycle 

( )j quantity associated with the /-th flight 

(j) angular coordinate for semi-elliptic surface crack, rad 


INTRODUCTION 


The NASA Dryden Flight Research Center (DFRC) B-52B launch aircraft has been used to cany various types 
of research vehicles by means of pylon hooks for high-altitude air-launch tests. The past air-launch vehicles 
include the following drop-test vehicles: the X-15 rocket plane (35,250 lb without drop tanks; 51,600 lb with drop 
tanks); FIL-10 lifting body (15,380 lb); highly maneuverable aircraft technology (FliMAT) vehicle (3,528 lb with 
4,000 lb adapter); the drone for aerodynamic and structural testing (DAST) vehicle (2,500 lb with 4,000 lb 
adapter); the solid rocket booster/drop test vehicle (SRB/DTV) (49,000 lb); Pegasus (Orbital Sciences Corp., 
Fairfax, Virginia) winged rocket (44,629 lb including 2,4001b adapter); and X-38 drop test vehicle (18,1001b). 
Each of the test vehicles is attached to the B-52B aircraft pylon through one L-shaped front hook and two identical 
real - hooks with the exception of the X-38. (The X-38 case used one front hook and one real' hook of a different 
pylon (ref. 1).) The L-shaped hook geometry will always induce tangential tensile stress concentration at the hook's 
inner curved boundary, which is the potential fatigue-crack initiation site. 

During the early stages of the SRB/DTV flight test program, the two old rear hooks (4340 steel) failed almost 
simultaneously during towing of the B-52B carrying the SRB/DTV on a relatively smooth taxiway (low-amplitude 
dynamic loading) after cancellation of the drop test due to unfavorable weather. Careful examination of the 
fracture surface of each failed old real' hook revealed that each hook had an existing microsurface crack at the 
critical stress points from where the microcrack rapidly propagated, resulting in the hook failure (ref. 2). The 
surface microcrack could have been initiated from the previous long period of flight-test stress cycling and surface 
coiTosion. Flad the hook failed during the takeoff run or during the captive flight, a catastrophic accident might 
have occurred. The potential for this type of accident underscores the need for reliable and accurate predictions of 
the fatigue life of pylon hooks. 

Using the half-cycle theory (ref. 3), Ko (refs. 4, 5) calculated the amount of crack growth at the critical stress 
point of each hook for each test flight using the actual random loading spectrum. Then, with certain assumptions, 
the Walker crack-growth equation was applied to formulate Ko first-order (ref. 5) and Ko second-order aging 
theories (refs. 6, 7) for predicting the number of remaining flights for the safe flight tests. Those theories predicted 
far fewer safe flights than those calculated from the classical aging theory. The Ko theories, however, were found 
to lose accuracy as the number of flights increases due to the neglected higher-order terms which become larger 
due to growing crack size. 

The safety of flight tests using aircraft pylon hooks to carry any drop-test vehicle hinges upon the structural 
integrity of the pylon hooks. It is, therefore, of vital importance to develop highly accurate aging theories to set the 
limit of safe flight-test life span for airborne critical structural components such as B-52B aircraft pylon hooks. 
This report presents two new aging theories for establishing the safe flight-test life span of B-52B hooks. The 
results are compared with the earlier Ko first- and Ko second-order aging theories, and also with the conventional 
aging theory for predicting the number of safe flights. 
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B-52B LAUNCH AIRCRAFT PYLON HOOKS 


Figure 1 shows the B-52B aircraft pylon carrying a particular store of Pegasus winged rocket through one front 
hook and two identical rear hooks before takeoff. Figure 2 shows the geometry of the front and rear hooks, with the 
location of the critical stress point for each hook indicated. Figures 3 and 4, respectively, show the tangential 
tensile stress distributions along the inner boundaries of the front and rear hooks calculated from the finite-element 
linear elasticity analysis (ref. 2). Note the tangential stress concentration at the critical stress point of each hook. 
Based on the finite-element stress analyses, the relationships between the hook loads { V A (front hook), V BL (left 
rear - hook), V BR (right rear hook) } (in lb) and the tangential tensile stresses { fi| , cr 7 , (7 3 , } (in ksi) at the respective 
hook critical stress points { 1, 2, 3} were previously established as (ref. 2) 

Front hook: 

<Tj = 7.3522 x 1(T 3 V a (1) 

Left rear - hook: 

a 2 = 5.8442 x l(T 3 V gL (2) 

Right rear hook: 

a 3 = 5.8442 x lCT 3 V g ^ (3) 

During the flight tests, the loading spectra of { V A , V BL , V BR } are obtained from the outputs of strain gages 
attached in the vicinity of the critical stress points of the hooks. The random cycling spectra of the critical point 
stress { <7 1 , (7 2 . fi 3 } are then calculated from equations (1) through (3) respectively. 

Figure 5 shows the random loading spectra for the three hooks during the takeoff and ascent of the B-52B 
carrying the SRB/DTV. Note that the magnitude of hook loading is more severe during the takeoff run than during 
airborne cruising. Figure 6 shows similar data during the landing of the B-52B carrying the SRB/DTV after an 
aborted flight test. Based on the data in figure 6, the peak hook loads at the moment of touchdown are listed in 
table 1 below. 


Table 1. Peak touchdown hook loads; SRB/DTV. 


Proof hook load, lb Touchdown hook load, lb Fraction of proof hook load 


- 36,520 

^A 

= 25,600 

/ = 

= 0.7 

= 44,110 

Vbl 

= 29,000 

/ = 

= 0.6574 

= 44,230 

V BR 

= 25,500 

/ = 

= 0.5765 
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From figures 5 and 6, it is seen that landing consumed far more fatigue life than the takeoff run. 

Figures 7-12 show portions of the random stress cycling spectra of { o^, a 2 - cs^} of the front and rear hooks 
carrying the SRB/DTV (figs. 7, 9, and 1 1), and the corresponding fatigue crack growth curves (figs. 8. 10, and 12). 
The fatigue crack growth curves shown in figures 8, 10, and 12 were calculated using the half-cycle theory to 
estimate the damage done by each half cycle of the random stress cycle spectra shown in figures 7, 9, and 1 1 
(refs. 3, 4), the only data available. Note that taxiing and takeoff consumed a great deal of hook service life because 
of higher crack growth rate due to high amplitude stress cycling because of ground surface effect. The crack 
growth rate then slowed during cruise because of the smooth air and the absence of ground surface effect. 


CONVENTIONAL AGING THEORY 

For each airborne structural component (e.g., B-52B pylon hook), let a 1 ’ be the initial fictitious crack size at 
the critical stress point associated with the static proof load test, and let a° c be the corresponding operational crack 

n q 

size associated with the peak operational load (a fraction of the proof load), then , a c } can be calculated 
respectively from the following crack tip equations based on fracture mechanics: 


P Q 
a' = - 

c 7t 


f K, 


\AM k g j 


0 Q 

a C = - 

c K 


f K, 


\AM K fa ) 


; /< i 


( 4 ) 


( 5 ) 


where K, is the critical stress intensity factor, a is the value of { < 7 ?,< 7 ~,<x! } at the respective critical stress 

c I 2 J 

point calculated respectively from equations ( 1 )— (3) using the respective hook proof load values 

P P P * 

{ V A , V BL , V BR } ; /( 7 (/■< 1) is the operational stress at the critical stress point; A is the crack location parameter 
(A = 1 for through-thickness cracks, A = 1.12 for both surface and edge cracks); and M K is the flaw magnification 
factor ( M K = 1 for very shallow surface cracks, M K = 1.6 when the crack depth approaches the back surface). 
Finally, Q is the surface flaw shape and plasticity factor for an elliptic surface crack (length 2c, depth a) (fig. 14, 
ref. 4) and is described as 



_ 

2 

( *\ 

Q = 

E(k) 

- 0.212 

(7 

l°r) 


where oy is the yield stress, and E(k) is the complete elliptic function of the second kind defined as 

u 

E W = f 2 a/ 1 -k 2 sin 2 0 
J o 


(6) 


( 7 ) 
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Because / < 1, the size of a° ( . is larger than that of a P (equations (4), (5)), the crack size differential 
(a° c - a P ) will be the limit of available crack growth for setting the number of safe flights. 

If Z\fl| is the amount of fictitious crack growth from the size a p during the first test flight, then the 
conventional age equation for estimating the available number of flights fj for a particular airborne critical 
structural component is given by 


F X 



Aa^ 


(8) 


Equation (8) assumes that the amount of crack growth per flight remains constant, and therefore, after each test 
flight, the number of safe flights will decrease linearly by just one flight. This assumption is highly inaccurate 
because the crack growth is progressive. As such, equation (8) is nonconservative and overpredicts the number of 
safe flights. Thus, for safe flight tests, more conservative and accurate aging theories are needed. 


CRACK GROWTH CALCULATIONS 


For constant-amplitude stress cycling, the fatigue crack growth is described by the following well-known 
Walker crack growth rate equation: 

^ - OTma x)'"( 1 - Rf= C(AK)"\ 1 -Rf~ m (9) 

where, C, m, n are material constants, is the maximum stress intensity factor, AK is the stress intensity 

amplitude, and R is the stress ratio given respectively by 


*max = AM k ° max 


na 


AK = AM k «J max -cr min ) 


na 


G 

_ nnn 


(10) 

( 11 ) 


(12) 


where G and a • are respectively the maximum and the minimum stresses of constant amplitude stress 

ITlaX 111 1 ll A J A 

cycles. 


For the random stress cycling, the half-cycle theory (refs. 3-7) may be used for the fatigue crack calculations. 
The half-cycle theory assumes that Walker crack growth equations (9) — (12) are valid to describe the damage done 
by each half cycle of the random loading spectrum. Under such an assumption, the amount of crack growth £>«■ 
caused by the i-th (i - 1 , 2, 3, ... ) half cycle of the random loading spectrum may be written from equation (9) by 
setting da = da . and dN = 1/2. Namely: 
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(13) 


5a = -C 
1 2 


,). 


"( 1 -*,■)"= 


with 


(^max)i 





(14) 


(^),.= AM, 






i - 1 


Q 


(15) 


R. = 






(16) 


where the subscript i (i = 1, 2, 3, . . . ) is associated with the i-th half cycle, and a- | is the crack size at the end of 
the (i - l)-th half cycle. 

If N, is the total number of random stress cycles induced by the /-th (/ = 1, 2, 3, . . . ) flight, then the amount of 
crack growth Aa, caused by the /-th flight may be calculated from 


Aa i = 



i = 1 


(17) 


In equation (17), the calculation of the right-hand side can be carried out through the use of special computer 
programs which sweep through all of the random stress-cycle data to pick up the values of {(< 7 max )„ (< 7 min ).} 
and calculate 5a i [equations (13)— (16)] for each half cycle, and then sum 5a - over 2 N, half cycles (not Nj) to 
obtain the value of Aa^ . 

It must be mentioned that the predictions of fatigue life by the half-cycle theory compare well with the 
experimental fatigue data (ref. 3). Figure 13 illustrates the method of summing up the crack growth increments 
caused by the random stress cycling when the half-cycle theory is used (refs. 5-7). 


EQUIVALENT-CONSTANT-AMPLITUDE AGING THEORIES 

The purpose of establishing the equivalent-constant-amplitude aging theories is to use the first-flight data to 
predict the number of future safe flights when the crack growth data for subsequent flights are not available. The 
approaches of the new aging theories are described below. 

First, the random loading spectrum obtained from the first-flight data is represented with an 
“equivalent-constant-amplitude loading spectrum”; then the Walker crack growth equation (9) is used to formulate 
the new aging theories for the calculations of the number of future safe flights. The equivalent-constant-amplitude 
loading spectrum is assumed to cause the same damage effect as the actual random loading spectrum based on the 
following assumptions: 
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1. The equivalent-constant-amplitude-loading spectrum induces the same amount of crack growth as 
Act i (/ = 1, 2, 3, ... ) caused by the random loading spectrum of each flight. 

2. The equivalent-constant-amplitude loading spectrum has the same number of stress cycles 
N, (/ = 1, 2, 3, . . . ) as the random loading spectrum of each flight. 

3. The values of { (7 max , R, N t } for the equivalent-constant-amplitude loading spectra remain the same for all 
flights. 

The equivalent-constant-amplitude aging theories will be formulated by the following two approaches: the 
“discrete” aging theory and the “closed-form” aging theory. 

Discrete Aging Theory 

For the discrete aging theory using an “equivalent-constant amplitude loading spectrum,” the Walker crack 
growth rate equation (9) is written in finite-difference form with each flight duration as a finite-difference interval. 
Then, the crack growth Aa, of the /-th flight is expressed in terms of the crack growth Au t of the first flight in the 
formulation of the discrete aging theory. 

1. Crack Growth 

Writing K x [equation (10)] for the (/- l)-th flight in the following form: 


K 


max 


AM k ° max 



( 18 ) 


where a, j is the crack size at the end of the (/ - l)-th flight. Then the Walker equation (9) may be written in the 
following finite-difference form to relate the number of stress cycles N, to the corresponding crack growth Aa, 
for the /-th flight (refs. 5-7) as 


( 


Acij = C 


AM ^ma xj 


(1 -R) Nj 


For the first flight, a { _ j = a j | = a’. , and equation (19) takes on the following form: 


(19) 


Aa^ = C 


AM k ° 'max 


\m 


m 

( \-R)\a P c ) 2 N l 


(20) 


Because the values of {<7 x , R, N ] } are assumed to remain the same for all flights, Zkq is proportional to 

m 

(a j | ) 2 only, and equation (19) written for the /-th flight may be divided by equation (20) written for the first 
flight to relate Acij of the /-th flight to Au | of the first flight. Namely, 
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Z\«l 

Aa^ 


i, 


= 1 


Aa 2 
Aa j 


m w 


'a | ^ 

2 

f Aa 1 

1 + 

P 


P 

V a c ) 


^ a c 2 


Aa^ 
Aa j 


v«F 


1 + 


ziflj + Zifl'2 


w m 


Aa 4 


2 

f Aa j + Aa 2 + Aa 2 

1 -i- 

Aa^ 

P 


p 

Ka c ) 


a c } 


( 21 ) 


( 22 ) 


(23) 


(24) 


Aa l 

Aa j 


m m 


(a ) 

a l- 1 

2 

Aa^ + Aa^ + Aa ^ + . . 

i 

sT 

+ 

p 

V J 



J 


(25) 


zlfl * 
F 

Aa j 


(a * ) 

2 

F -1 

_ 

p 


V a c ) 



1 + 


1 + "t" + - - - + 21^7 ^ + - - 


+ Aa 


m 


F -1 


(26) 


It is seen that the crack growth Aa j for the Z-th flight is expressed in terms of the crack growth of all previous 
flights (Aa^, Aa 2 , Aa-j, Aa 4 , . . . , Aaj _ ] ) , and that the amount of crack growth is progressive, meaning the 
amount of crack growth Aa t will continue to increase with the increasing number of flights 
(i.e., Aa^ < Aa 2 < Aa 3 < . . . < Aaj < . . . ). 

In light of equations (21)-(26), if the crack growth Aa | of the first flight is known (i.e., calculated from 
equation (17) from the first-flight data), then the crack growths (Aa 2 , Aa 2 , Aa 4 , . . . , Aaj, . . . , Aa ^.*) for all 
subsequent flights can be calculated. This is the essence of the discrete aging theory. 


For mathematical simplification, the right-hand sides of equations (21)-(26) have been expanded to the first- 
and second-order terms. The Ko first-order (ref. 5) and Ko second-order (refs. 6, 7) explicit aging equations were 
developed for the calculations of the number of safe flights {Fj , F i } . Those explicit aging equations based on the 
Ko first- and Ko second-order aging theories are taken from references 5-7 and are referred to in the appendix. 
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2. Number of Safe Flights 


When the crack growth summation (Aa^ + Aa 2 + Aa 2 + ... + Aa ^ + Aa^,*) reaches the size of the 
available crack differential ( a° c - a p ) , one can write 

Aa^ + Aa 2 + Aa 3 + . . . + Aa * ^ + Aa^- t . = a° c - a P = F jX Aci^ (27) 

* 

|< F terms >| 


Equation (27) is an implicit aging equation based on the discrete aging theory for finding the number of safe 

* 

flights F . By substituting equations (21)-(26) into equation (27) with Aa t of the first flight calculated from 

* 

equation (17). the number of safe flights F can be calculated. 

* 

In the computational process to find the value of F from equation (27), the number of terms of summation of 
Aa, on the left-hand side was successively increased step-by-step until the value of summation of the left-hand 

O D 

side approached the value of {a c - ) on the right-hand side. 

Closed-Form Aging Theory 

For the closed-form aging theory based on the “equivalent-constant-amplitude loading spectrum”, the aging 
equation will be formulated through the integration of the Walker crack-growth-rate equation (9). The closed-form 
aging theory will then be used to check the accuracy of the discrete aging theory. 

1. Crack Growth 

To formulate the closed-form aging theory, equation (10) is substituted into the Walker crack-growth-rate 
equation (9) to yield the following form: 


da 

dN 


AM K ° 'max 



(1 -/?)"(«) 


m 

2 


(28) 


which, after rearranging, becomes 


da = C 


AM K a max 


n 


Q 


(1 -R) dN 


(29) 


Because the values of { c max , R} are assumed to remain unchanged for all flights, equation (29) may be easily 
integrated from a p to arbitrary crack size a for the left-hand side, and from 0 to the corresponding number of stress 
cycles N for the right-hand side as 


10 



(30) 




= C 


m 






(1 -R)N 


which may be re-written in alternative form as 


a 



C 


m 
2 ' 


AM K G max 




2 - m 


(i 


R) n N 


(31) 


Equation (31) is to be used to plot the crack size a as a function of number of stress cycles N. 

Also, integration of equation (29) over each flight interval with N, number of stress cycles, yields the 
following “powered” crack-growth equation for each flight: 

First flight: 





AM K G max 


n 


Q 


(1 -R) N j 


Second flight: 




m 
2 - 


= c 


m 


-1 


AM K G max 




(1 -R) N 2 


(32) 


(33) 


/-th flight 




= C 


m 


-1 


AM K °. max 


\m 


(1 -R) N 


(34) 


Fj -th flight (Final flight) 


(1 -R) n N* ; a * = a° (35) 

F 1 F j C 

The cai licr assumption that the loading spectra of all flights have the same number of stress cycles: 


a * 

F , - 


Cl * 

F , 


= c 


m 
2 ' 


AM K ° max 
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( 36 ) 


Nl = N 2 = N 3 = ... =N,= ... =N p * 
gives the following “powered” crack-growth relationships between all flights: 



2. Number of Safe Flights 

* Q 

By adding all the crack-growth equations (32)-(35) from flight / = 1 to / = F , , and noting that a * = a , 

i p c 

there results: 1 


4V 1 -k 


m 
' 2 


(38) 


= C I - - 1 


AM v o 

K max, 


\n 


(1 -R) [N l + N 2 + . 


. +N,+ 


. +N 


Dividing equation (38) by equation (32) results in an explicit aging equation for the calculation of the number 
of safe flights F ^ : 


|< F terrns- 


A^ 1 +A^ 2 + A^3 + ...+A^ / +... +A^ . 

F i 


N, 


1 _m j 

4) 2 -f«^ 2 


m 


(39) 


a P c 1 2 - I a, 


1 — 


which, in light of the assumed condition (36), becomes 


f l^l 

*1 



(40) 


* 

which is the aging equation for calculating the number of safe flights F j based on the closed-form aging theory. 
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Alternatively, writing equation (30) for the flight age limit (i.e., a = a° c , N = N toJa j), we have 


• m | _ m 

KJ ' 5 

* 

Because of the assumption N. ta/ = F^xN^ 
equation identical to equation (40). 

It must be mentioned that for the value of crack size a j = a p + Aa^{l = 1 ) for the first flight, Aa j is to be 
calculated from equation (17) using the half-cycle theory as mentioned earlier. 

EQUIVALENT STRESS AMPLITUDES 

In the previous formulations of the equivalent-constant-amplitude aging theories, { (7 max , f?, A ; } were 
eliminated through division processes and, therefore did not appear in the final aging equations. Thus, there is no 
need to know the actual magnitudes of { <7 , R, N A . 

In fact, each random loading spectrum could be represented by a multiple number of equivalent- 
constant-amplitude loading spectra based on the selections of the values of { <7 , R, N,} . 

For the purpose of visualization, we will calculate the values of { cr , <7 min , R } , and plot against the random 
loading spectra flight data for comparisons. 

To determine the maximum stress <7 max and the mean stress <7, of the equivalent-constant-amplitude loading 
spectrum, equation (30) may be rewritten in the following form: 


AM K° 'max, 




(1 -R) N 


total 


(41) 


, division of equation (41) by equation (32) will yield the aging 


(I"*) = 


C(m - 2) 


1 \Q 

{ AM k aJ k j 


, m , m 

P f~2 ( 

\m \a c \ z - \ fl, 


A, 


(42) 


For the equivalent-constant-amplitude load spectrum, the following relationship holds: 


R = 


%ax 

(43) 

2(7 


s 1 

(44) 

%ax 



Substitution of relationships (43) and (44) into the left-hand side of equation (42) yields: 
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Change No. 1. April 1, 2004 


«W)‘ 


l — — — 


a 

V max/ 


( 2) 1 ” 
C(m - 2) 


[ AM k 




(45) 


In equation (45), the right-hand side is known, and the left-hand side has two unknowns: <7 max and o s . Thus, 
there are multiple combinations of c max and o s for finding the equivalent-constant-amplitude load spectrum. The 
most practical way is to choose the value of the mean stress a s of the random loading spectrum by inspection, and 
then calculate the equivalent maximum stress <y from equation (45). 


NUMERICAL EXAMPLE 


The numerical example is to demonstrate how to calculate the number of safe flights for a flight test program. 
The example chosen is the NASA B-52B carrier aircraft pylon hooks (one front hook, two rear hooks) (figs.l, 2) 
carrying the Pegasus winged three-stage solid rocket (44,629 lb) up to high altitude (approximately 40,000 ft) for 
air-launching, firing and sending the payload into orbit. 

Material Properties 

The material properties of the B-52B pylon front and real - hooks are listed in table 2. 


Table 2. Material properties of B-52B pylon hooks. 


Part name 

Material 

°U 

ksi 

<*Y 

ksi 

T U 

ksi 

Kr 

L c 

ksi Vim 

C 

in. ,, ■ r . — 
(ksiVm.) 

cycle 

m 

n 

Front hook 

Inconel 718 

175 

145 

135 

125 

9.220 X 10“ 12 

3.60 

2.16 

Rear hooks 

AMAX MP35N 1 " 

250 

235 

141 

124 

2.944 X 10 -11 

3.24 

1.69 


* Inco Alloys International, Inc., Huntington, West Virginia, 
t H. C. Starck, Inc., Cleveland, Ohio. 


Input Crack Data 

The fictitious surface crack at the critical stress point of each hook is assumed to be semicircular in shape 
(c = a, fig. 14). This is based on the observations of the past failed B-52B pylon two rear hooks. The crack 
propagation initiation sites of the failed hooks were almost semi-circular surface cracks of depths a = 0.03 1 in. and 
a = 0.038 in. respectively for the inboard and outboard hooks (ref. 2). Note from figure 14 that the semicircular 
surface crack has the relatively high value of Q. Other crack parameters used were: A = 1.12 for the surface crack, 
and M K = 1.0 for the high (hook depth) to (crack depth) ratios. 
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Before the test flights, the approximate number of safe flights must be estimated. Because the actual 


random loading spectrum for each hook is not available for the calculation of the amount of crack growth Aa^ for 

Act , 

the first flight, the value of — - for each hook may be based on the earlier test data of SRB/DTV because of the 


weight proximity of the two vehicles. The SRB/DTV (49,000 lb) is only 8.92 percent heavier than the Pegasus 


winged rocket (44,692 lb). 


Aa, 


P O AALl i 

The initial crack size a c , the operational crack size a r , and the crack growth used in the aging analysis 

a c 

for the case of the Pegasus winged rocket are listed in table 3. 


Table 3. Proof and operational crack sizes for B-52B pylon hooks carrying Pegasus 
winged rocket (44,692 lb). 


Hook name 

Hook proof 
load, lb 

P ■ 

a c , in. 
(/= 1-0) 

O . 

a c m. 
(/= 0.6) 

Auq 

a c 

Front, V A 

36,500 

0.1247 

0.3465 

0.01814 

Rear left, V BL 

57,819 

0.0774 

0.2151 

0.00761 

Rear left, V BR 

57,819 

0.0774 

0.2151 

0.00989 


In table 3, the value of initial crack size a P and the operational crack size a° c are calculated respectively from 
equations (4) and (5). The value of a° c was calculated based on the assumption that operational peak stress is 
60 percent of the proof stress (i.e.,/= 0.6). This/value is slightly higher (more conservative) than the/values in 
the SRB/DTV case, for which the / value for the front hook is /= 0.5450, for the left real - hook,/ = 0.5946, and 
for the rear right hook,/ = 0.5986 (refs. 4, 5). 

Aa^ Aa^ 

Note that in table 3, the crack growth of the right rear hook is 1.3 times larger than of the rear left 

P P 

a c a c 

hook. This is based on SRB/DTV flight data which show that the crack growth rate for the right rear hook is 
1.3 times larger than that of the left rear hook as a result of the combined aerodynamic, inertial, and static loads. 

Equivalent-Constant-Stress Amplitudes 

The typical duration of a B-52B air-launch test flight is about 50 minutes (in the SRB/DTV case), 
and the random loading spectrum has 4 cycles per second. Therefore, the number of stress cycles is 
Vj = 12,000 (= 50 X 60 X 4) cycles per flight. This value will be used in generating the equivalent- 
constant-amplitude loading spectrum. Table 4 shows the calculated value of {a m ax ,c min , R} associated with 
the equivalent-constant-amplitude load spectrum for given N j , Aa j / a P c and cr y . 
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Table 4. Stress values for the equivalent-constant-amplitude loading spectra; N j = 12,000 cycles: Aa^/a p and 
a s given. 


Flook 

Aa j 

— (Given) 

fl c 

o s , ksi (Given) 

%ax ’ ksi (/) 

°min- ksi 

R 

Front ( V A ) 

0.01814 

108.000 

115.406 (0.43) 

100.594 

0.8717 

Rear left ( V Hl ) 

0.00761 

127.000 

130.157 (0.50) 

123.843 

0.9515 

Rear right ( V BR ) 

0.00989 

132.500 

136.051 (0.53) 

128.949 

0.9478 


Figures 15-17 compare the actual random loading spectra (using available data of SRB/DTV) and the 
associated equivalent-constant-stress amplitude stress cycles for the three hooks based on the data listed in table 3. 
Notice that the / values for the equivalent-constant-stress amplitude stress cycles to cause the same amount of 
damage as the random loading spectra are only in the range of /= 0.43-0.53 (much less than/= 0.6). 

Number of Safe Flights 

Figures 18-20 show normalized crack sizes (Aa^/Aa^) plotted as functions of the number of the flight (/) for 
the front and the real - hooks based on the different aging theories. The terminal point of each crack growth curve 
gives the predicted number of flights. Notice that the conventional aging theory gives horizontal lines (figs. 18-20) 
because of the assumption, Aa l = Aa^ . As the aging theory is improved, the slopes of the crack growth curves 
become steeper and the predicted number of flights decreases. The crack growth curves based on the discrete aging 
theory and the closed-form aging theory are practically coincidental. The two theories predicted identical and the 
lowest number of flights. 

Figure 21 shows the crack growth curves in the a - N space for the front and rear hooks based on the 
closed-form aging theory [equation (31)]. For each hook, the intersection of the crack growth curve and the 
corresponding operational crack limit line ( a = a° horizontal line) will give the total number of available stress 
cycles. After dividing the total number of cycles by 12,000 cycles per flight, one obtains the number of remaining 
flights Fj = 39, Fj = 100, Fj = 77 flights respectively for the front hook, rear left hook, and rear right hook. 

Figures 22-24, respectively, show the crack growth curves for the front and the two rear hooks in the 
a - flights space. The crack growth curve based on the classical aging theory for each hook is a straight line 
because of the assumption that the amount of crack growth Aa, for every flight remains the same. The crack 
growth curves based on the other aging theories are nonlinear. The intersecting points of the crack growth curves 
and the associated horizontal operational crack limit line ( a = a" horizontal lines) will give the number of 
remaining flights {Fj, F^, Fj, F* F* } based on the different aging theories. The values of 

— ~ 4 : 4 = 

{Fj, Fj.Fj, F,Fj } are indicated in the figures. 
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Table 5 summarizes the number of remaining flights {F |; F|, F|, F* F* } predicted from different aging 
theories for the case of B-52B hooks carrying the Pegasus launching vehicle. 


Table 5. Number of safe flights predicted from different aging theories for the case of B-52B hooks carrying the 
Pegasus winged rocket (44,692 lb). 


Hook 

F l 

(f x /f\) 

F\ 

{F x /f\) 

F 1 

Cf x /f\) 

* 

F 

* 

F 1 

Front ( V A ) 

98 

(2.51) 

53 

(1.36) 

45 

(1.15) 

39 

39 

Rear left ( V BL ) 

234 

(2.34) 

130 

(1.30) 

110 

(1.10) 

100 

100 

Rear right (V BR ) 

180 

(2.34) 

100 

(1.30) 

85 

(1.10) 

77 

77 


In table 5, the predictions from different aging theories are also normalized by the closed-form solution Fj 
(inside the parentheses). Note that the classical aging theory predicts the number of remaining flights by more than 
twice that predicted by the closed-form (or finite-difference) aging theory. 

CONCLUDING REMARKS 


New aging theories have been developed to accurately estimate the number of safe flights for airborne 
structural components. The highlights are listed below. 

1. The newly-developed aging theories take into account the progressive crack growth problem, and, 
therefore, predict far fewer safe flights of airborne critical structural components than that predicted from 
the conventional aging theory due to the inclusion of all the higher-order terms. 

2. The new aging theories and the previous Ko first- and Ko second-order aging theories agree well at a low 
number of flights. However, as the number of flights increases, the two previous aging theories slightly 
oveipredict the number of remaining flights in relation to the two newly-developed aging theories which 
include all the higher-order terms. 

3. Both the discrete and the closed-form aging theories predict a practically identical number of safe flights. 

4. The closed-form aging theory provides an explicit expression of the aging equation for the calculations of 
the number of safe flights. 

5. The discrete aging theory allows step-by-step updating of the number of remaining safe flights for 
geometrical visualization of how the number of remaining safe flights are consumed. 

6. The front hook has the shortest flight-test life span, and therefore is the most critical structural component 
in determining the limit number of flights for safe flight tests. 

7. In the case of the Pegasus winged rocket, the number of available safe flights for the flight tests is 
39 flights. 
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APPENDIX 


KO FIRST-ORDER AND KO SECOND-ORDER AGING THEORIES 

This appendix is based on the results presented in references 5, 6, and 7. 

Under the assumption that {<7 m , R, NA remain the same for all flights, and while the value of 


A a i+Au^ + Ar/^ + ... 


is still small (i.e., 


Aa y + Aci^ ■+- Aaj + . . . 


P p 

a c Aa, a c 


« 1) during the earlier stage of flight tests 


C LAU 1 C 

(i.e., / is not large), the crack growth for each flight [equations (2 1)— (26)] may be expanded up to 

Aa j 

second-order terms as follows (refs. 5, 6, 7): 


(A-l) 


Aa,) 2 ( m 

1 + — 1 = 1 + 1 - 


m^ a l j 1 m A a \ + \^(\ _L^| m ^ a l 2 + 

2 jp 2 2 p V m) 2 p 


(A-2) 


a 2 Y 


Aflj + Aa-, |2 


t,f Aa-,31 2 

- 1 + — 


(A-3) 


- 1+2 m ^ a \ y 1 m ^ a l + 2 2 fl 1 ” + 


2 of A A 2 of 


'"A 2 a p c 


a 3 | 2 


Auq + Aa-, + AaY 2 


Aa ,( Aa ~ Aa-, 

1 + 1 + — = + — 2 

P I Aflj Aciy 

« c v 


(A-4) 


- 1+3 m y 1 m ^ a i +3 2 fi 13 m ” + 

+ 2 p 22p + vm72p + 

v a c 7L v a c 7J v \ a c J 
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(A-5) 


m m 


Ac^ = 

i a l-\ 


2 

1 Aci | + . . 

■ +Aa i-\ 


P 

^ a c 

9 



y 


Aa , f Aa~. Act-, 

1 + — -| — - + — l + . . 

Act, Aa , 


Aa 


/- 1 




'1 


Aa 


1 )_ 


= ! + (/-!) 


m Aa l 

"l-I 

2 

m Aa l 

i 

A 

i 

+ 

m Aa l 

2 p 

2 p 

2 p 

V a c J 


V a c J. 

j 

K a c ' 


Ko First-Order Aging Theory 

Aa, 

The first-order expansion of [equations (A-l)-(A-5)] gives the following first-order aging formula 

Aflj 

for the calculation of the remaining flights F | (ref. 5). 


Fl 



mAa j 



mAa j 



^t 


\ 

-1 

/ 


(A-6) 


Notice that the number of remaining flights F \ based on the Ko first-order aging theory is expressed explicitly 
in terms of Fj , the number of remaining flights based on the classical aging theory. 


Ko Second-Order Aging Theory 

By retaining up to the second-order terms of expansions in equations (A-l)-(A-5), one obtains the 
following Ko second-order aging equation for the calculations of remaining flights F \ expressed explicitly in term 
of Fj as follows (refs. 6,1). 


Fj = B + D -- 
1 3 

where both B and D are functions of Fj and are given by 


5 = 3 j- £ + 1 

l£ V 

A| 2 V 

4 

27 

F> = 3 L@_ 

o 

V 2 aJ 

4 

27 


(A-7) 


(A-8) 


(A-9) 


where 
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(A- 10) 


1 2 , 
a = ^Oq-p ) 


P = ±(2p 3 -9pq + 21r) 


(A- 11) 


where p, q, r arc the coefficients of the following cubic equation 


~ 3 ~ 2 

F\+ pF \+ qF j + r = 0 


(A- 12) 


and are given below. 


3 , 3 , 0 a c 

p = 1 — m + Z 

2 (m —1)1 2 Zlflj 


(A- 13) 


q = I 1+-L_ ™-2 + i 

^ ... 1 ^ A .. A .. 


m - 1 2 Aa^ m\Aa^ 


(A- 14) 


m(m - l)\Aa ^ 


(A- 15) 


Notice that Fj (remaining safe flights based on the conventional aging theory) appears only in the 
expression for r in equation (A- 15). 


At relatively low number of flights (i.e., / is not too large), the term 
Aci , + An-, + do, + . . . I 


A« j + + Au^ + . . . 


is still small 


i.e., «1 , the Ko first- and Ko second-order aging theories could give reasonably 

l a P c J 

accurate number of remaining flights. However, the two theories start to lose accuracy as the number of flight (/) 


increases. 
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Figure 3. Distribution of tangential stress c r along inner boundary of front hook; V A = 10, 000 lb. 



Figure 4. Distribution of tangential stress a r along inner boundary of rear hook; V BL - V BR - 17, 179.53 lb. 
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Figure 5. Random loading spectra for B-52B hooks carrying SRB/DTV (captive flight); takeoff run and airborne. 
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Figure 6. Random loading spectra for B-52B hooks camying SRB/DTV (captive flight); landing and taxiing. 
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Figure 7. Typical random loading spectrum for front hook ( V A ) critical stress point 1 (cij); B-52B carrying 
SRB/DTV. 
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Figure 8. Crack growth curve for front hook ( V A ) critical stress point 1 ( ); B-52B carrying SRB/DTV. 
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Figure 9. Typical random loading spectrum for left rear hook (V BL ) critical stress point 2 (c 2 ); B-52B carrying 
SRB/DTV. 
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Figure 10. Crack growth curve for left rear hook (V BL ) critical stress point 2 (cr 9 ); B-52B carrying SRB/DTV. 
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Figure 11. Typical random loading spectrum for right rear hook ( V RR ) critical stress point 3 ( (J-, ); B-52B carrying 
SRB/DTV. 
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Figure 12. Crack growth curve for right rear hook ( V BR ) critical stress point 3 ( <r 3 ); B-52B carrying SRB/DTV. 
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Figure 13. Graphic evaluation of crack growth increments for random loading spectrum using half-cycle theory. 









Figure 16. Representation of random loading spectrum with equivalent-constant-amplitude loading spectrum; 
rear left hook ( V BL ). 



Figure 17. Representation of random loading spectrum with equivalent-constant-amplitude loading spectrum; 
real - right hook ( V BR ). 
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Figure 18. Crack growth ratio (A a, / Aa j ) curves calculated from different aging theories; front hook ( V A ). 



Figure 19. Crack growth ratio (Aa, / Aa | ) curves calculated from different aging theories; rear left hook ( V BL ). 
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Figure 20. Crack growth ratio (Aa, / Aotj ) curves calculated from different aging theories; rear right hook (V BR ). 
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Figure 2 1 . Crack growth curves in a - N space calculated from closed-form aging theory for the front hook ( V A ) 
and two real - hooks {V BL , V BR ). 
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Figure 22. Crack growth curves in a - flights space calculated from different aging theories; front hook. 



Figure 23. Crack growth curves in a - flights space calculated from different aging theories; real - left hook. 




Figure 24. Crack growth curves in a - flights space calculated from different aging theories; rear right hook. 
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